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ABSTRACT 

Pulsar Timing Arrays (PTAs) use high accuracy timing of a cohection of low timing 
noise pulsars to search for gravitational waves in the microhertz to nanohertz frequency 
band. The sensitivity of such a PTA depends on (a) the direction of the gravitational 
wave source, (b) the timing accuracy of the pulsars in the array and (c) how the available 
observing time is allocated among those pulsars. Here, we present a simple way to cal- 
culate the sensitivity of the PTA as a function of direction of a single GW source, based 
only on the location and root-mean-square residual of the pulsars in the array. We use 
this calculation to suggest future strategies for the current North American Nanohertz 
Observatory for Gravitational Waves (NANOGrav) PTA in its goal of detecting single 
GW sources. We also investigate the affects of an additional pulsar on the array sensi- 
tivity, with the goal of suggesting where PTA pulsar searches might be best directed. 
We demonstrate that, in the case of single GW sources, if we are interested in maximiz- 
ing the volume of space to which PTAs are sensitive, there exists a slight advantage to 
finding a new pulsar near where the array is already most sensitive. Further, the study 
suggests that more observing time should be dedicated to the already low noise pulsars 
in order to have the greatest positive effect on the PTA sensitivity. We have made a 



web-based sensitivity mapping tool available at jhttp : //gwastro . psu . edu/ptasm 



Introduction 



Pulsar Timing Arrays (PTAs) (JFoster &: Backeiill99Cll ) are a practical means of observing gravi 
tational waves (GWs) associated with supermassive (> IO^Mq) b lack holes, or more generally GWs 



in the nanohertz to microhertz regime arising from any source (IJaffe &: Backer 



2006a 



2003 



Jenet et al. 



Sesana Sz Vecchid l2010al : lOlmez. Mandic. &: Siemens! I2OIOI ) . A PTA is a collection of low 



timing noise pulsars timed to high accuracy. In such a timing array, GWs may signal their presence 
through correlated disturbances in measured pulse arriva l times. Over the last dec ade, assessments 
of expected GW source strengths have remained steady (jSesana Sz Vecchid l2010al ) while the num- 
ber of pulsars that can be timed to high precision has increased, and instrumentation, observing 
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technique, and timing precision of known pulsars have all improved (jHobbs et al.ll201Cll ). Detection 
of GWs via pulsar timing is now within striking distance. Properly made strategic choices for 
observing and/or improving the timing precision of current array pulsars, and searching for new 
pulsars, can hasten the arrival of the day when the first detection is made and the field of GW 
astronomy is properly inaugurated. 

Strategic optimization depends upon your strategy. Do you want to detect single sources, or a 
stochastic background? Do you want to maximize the number of sources detected, or maximize the 
signal-to-noise ratio (SNR) of the sources that you do detect? Here we investigate the sensitivity of a 
PTA to a single source of GWs and investigate optimization strategies that maximize the number of 
sources detected. We describe the results of a preliminary investigation into strategies for deploying 
limited observing time among current array pulsars, and for searching for new pulsars to augment 



the current array. O ur results apply to the search for single sources of GWs(iFinn &: Lommen 



Yardlev et al.ll20ldl. not t o stochastic background searches (jJenet et al.ii2006b : iAnholm et al 
van Haasteren et al.ll2009l ). 



201C 



2008 



After an introduction of the PTA single-source sensitivity calculation and NANOGrav charac- 
teristics in §2, we discuss in §3 the consequences of adding pulsars to the array. Further, discussed 
in §4, is the preliminary study of pulsar timing allocation time. Finally, in §5, we discuss the 
implications that rise from this study and suggest future research. 



2. Calculating PTA Sensitivity 

Anticipated GW so urces in the nHz-//Hz band include nearby supermassive black hole b i- 



nary and triplet systems (jJaffe fc Backer 



2003 



Sesana &: Vecchio 



2010a 



Amaro-Seoane et al.ll20ld ). 



More speculativ e sources include a primordial stochastic background and bursts from cosmic string 



cusps and kinks ( Damour fc Vilenkinll200ll . l2005l : ISiemens et al.ll2006l : ISiemens. Mandic. Sz Creighton 



2007 



Qlmez. Mandic. fc SiemensI l20ld ). All these sources are expected to be isotropically dis- 



tributed on the sky. Correspondingly we choose to estimate the sensitivity of PTAs as the spatial 
volume within which a fiducial source would give an SNR greater than a fixed threshold. Such a 
measure can be made independent of the threshold (but not the fiducial source) by referring the 
calculated volume to a reference volume calculated for a reference array. In this way we define the 
relative overall sensitivity Voverall of a PTA as: 



^overall 



Etidik) 



(1) 



R 



where d is the distance out to which the PTA is sensitive to a source propagating in direction ki 
relative to the reference PTA denoted with a subscript R. The limit n is the number of directions 
along which sensitivity is being measured and can be chosen by the user according to the desired 



resolution of the resulting sensitivity map. In the maps we constructed in this paper we used 
HEALPi^o to construct pixels representing equal surface area, and n = 3072. 

The distance to a fiducial GW source in a fixed direction is inversely proportional to its 
amplitude SNR p as observed in a PTA. Correspondingly, we may use the amplitude SNR of a 
fiducial source (assuming we have measured the optimal SNR) as a surrogate for the distance of 
the source 



d{k,p) = d{k,po) — , 
P 



(2) 



where d{k,p) is the distance to the source when its anticipated SNR is p, its propagation direction 
is k, and po is the amplitude SNR of the fiducial source at distance d{k,po). 

Timing noise for typical PTA pulsars is typically white on timescales les s than 5-10 years 
and red on longer timescales (jJenet et al.l l2006al : iHobbs. Lvne. fc Kramerl |2006| ) . The white con- 
tribution to the timing noise is characterized by its root mean square (RMS) residual aj. For our 
approximate analysis here we ignore the red timing noise contribution and assume, for each PTA 
pulsar, white timing noise characterized by aj for pulsar j, and that a single aj characterizes the 
entire observation. 



With this timing noise approximation the contribution from pulsar j to the PTA power SNR 
in direction k is 



p^(/c,t) 






(3) 



where Tcwik, t)j is the anticipated GW contribution to the residuals in pulsar j from a GW source 
propagating in direction k at time t. The anticipated power SNR p^ of the whole array to sources 
propagating in direction k is the sum over the pulsars 



p\k,t) 



V P GW 



(k,t)j 



(4) 



where Up is the number of pulsars. Armed with this equation, we set off to find TGwik,t)j. 



The timing residuals tqw associated with TT-gauge GW metric perturbation /i+e_)_ + /ixe> 
may be written as 



TGwik,t)j = F+{k,nj)g+{t,Lj,k-hj) + F^{k,hj)gy,{t,Lj,k ■ fij), 



(5) 



where hj is a unit vector pointing to pulsar j, Lj is the distance to that pulsar, and k is the direction 
of propagation of the GW. F^' ^ are given by 



,nj) 



F+{k 



F''{k,hj 



^nje'^{k)hj, 
inje^(fe)nj. 



(6) 

(7) 



where h denotes the transpose of h, and e'+^ and e'^-* are the two independent gravitational wave 
polarization basis tensors, 



Im J+) _ Jm (x) 
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Functions 5+ and gx are integrals of /i+ and /ix as follows (JFinn fc Lommenl 12010 ): 

9{+/x)it, Lj, kj ■fv') = j /i+/x [t-{l + k- nj){Lj - A)j d\. 



(8) 
(9) 

(10) 



Note that we are using geometrized units where c = G = 1. We have essentially broken up 
' Tnw{k)j into terrns tha t depend on geometry (F's) and terms that depend on time (5's). Following 
([Finn &: Lommenl 1 20 lOl ) we assume that a function / exists for which 



df+/x{u)/du = h+jxiu). 
We can then do the integral as follows: 

/+/x(t) /+/x(t-(l + fc-n,-)L,) 



g^+/x){t,Lj,kj -h^) 



1 + k ■ hj 



1 -\- k ■ hj 



(11) 



(12) 



The first term is the so-called 'earth term', and the second the 'pulsar term' ljjenet et al.1 12004 ) . 
The pulsar term is delayed from the earth term by (1 + /c • nj)Lj which amounts to hundreds to 
thousands of years in most cases. For the moment, assume that we are dealing with burst sources 
whose length is shorter than our observation time (years) for which we only observe the earth term, 
and that we can thereby ignore the pulsar term. Later we will show that the result we derive here 
holds for continuous sources, when the pulsar term must be included, as well. 

TGwik,t,t)j depends on time, so what we desire is the time averaged SNR, p^, i.e. the time 
average of equation [H We square equation [5] and average over time and all possible polarizations /i+ 
and hx- The cross-term F~^{k,hj)F^(k,hj) vanishes when we average over all polarizations. Also 
by averaging over all polarizations we find g'^ = g'^. Finally, exploiting time translation symmetry 
we can write 



A^ 



F(+)2(fc,nj)+F(x)2(fc,nj) 
a2(l+fc.fi,)2 



(13) 



where A^ is a constant independent of the pulsar line-of-sight hj or the propagation direction of 
the GWs k. 



Combining equations [21 HI [6l [7] and [13] we have our principal result 
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(15) 
(16) 



The steps between equations 14 and 15 can be done for any k, but the resuh in equation 15 can 
be seen more readily by assuming the GW is travehng in the z direction which gives 
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and assuming an arbitrary pulsar direction fij = [nx,ny, Uz]. The numerator in the sum in equation 
14 becomes (after some algebra and trigonometry) (1 — n^)^/4. For an arbitrary GW direction, k, 
the numerator generalizes to {1 — {k ■ nj)-^)^/4. In other words, the quantity that matters is the 
projection of the pulsar direction vector onto the direction of propagation of the GW. 

Equation 16 will allow us to compare various PTAs independently of the details of the input 
source {g+ and gx as shown in equation [5]) and then putting this into equation [D 



-^^3/2 



^overall 



.iPHki) 



SLipK^O 



^^372'' 



(17) 



where p^ refers to the reference PTA, and, as in equation[Tl the n is the resolution of the calculation, 
i.e. the number of pixels in the map of p'^{ki). ^overall is what we are calling the 'volume sensitivity' 
and represents the ratio of the volumes to which two different PTAs are sensitive. 

A related quantity which we will utilize later, y{k), represents the comparison of the sensitivity 
of two arrays as a function of the GW propagation direction k, 



uCk) 



-T^3/2 



pHk) 



-^^3/2 • 



Pls(fc) 



(18) 



Note that rather than plotting i^{k) as a function of GW propagation direction k we will plot P(s) 
where s is the direction of the GW source in the sky, s = —k, and i>(s) = u{—k). 



The quantity 



P\k) 



pHk) 



maxr p2 (k) 



(19) 



lends itself to interpretation as the PTA antenna pattern: i.e., it is proportional to the signal power 
absorbed by the detector from a source propagating in the direction k, measured relative to the 
source direction for which the greatest power is absorbed. As with i'{k), P^(s) = P'^(—k). Figure 
1 shows "source- averaged" antenna pattern P^{s) for the NANOGrav PTA, whose member pulsars 
and their RMS timing residuals at the time of writing are provided in Table 1. We have also made 



available a web-based tool for computing sensitivity maps for an arbitrary array of pulsars at this 
URii. 




Fig. 1. — Antenna pattern -P^(s) (defined in the text following equation llOh of the NANOGrav PTA 
is shown in equatorial coordinates. In other words, this shows the sensitivity of the NANOGrav 
PTA to single GW sources as a function of location of the source. The diamonds show the location 
of the pulsars in the array. 



As evident from Figured] the sensitivity of the NANOGrav PTA is heavily biased to the region 
between 12 and 24 hours right ascension, yet fairly symmetric about the equatorial plane. Not 
surprisingly this region represents both the highest concentration of pulsars and also the location 
of the lowest noise pulsars. 

Rather than average over all possible GWs propagating in direction k we can restrict attention 
to any particularly interesting class of sources. Consider, for example, the radiation from a circular 
binary consisting of two supermassive black holes. For syster us like these the radiation is approxi- 
mately periodic over any reasonable observational timescale (jJenet et al.ll2006bl : ISesana &: Vecchio 
2ninbh and 



TGwik)j oc cosujt — cosuj t — {1 — k ■ nj)Lj 



(20) 



where uj is the GW angular f requency. Even in th e best of circumstances pulsar distances are 
known to no better than 10% (ICordes &: Lazidl2002l ). in which case the phase loLj is uncertain by 
many times 27r for oj of interest. If we average TGw{k)j over the typical uncertainty in distance 
we find the contribution owing to the term involving Lj, the pulsar term, vanishes and we are left 
with cos^wt which, averaged over time, tends to 0.5, a constant. As long as the uncertainty in 
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pulsar distance is greater than the light travel time over the duration of the observation these same 
considerations will hold for any source that a PTA can detect: i.e., averaged over the uncertainty in 
pulsar distances the pulsar term contribution to r^ will vanish and we can again ignore the second 
term on the right-hand side of equation [12] as we did to obtain the result shown in equations 16. 
Note that in the case of k ■ h = 1 the earth and the pulsar terms exactly cancel and the sensitivity 
p'^{k) = as is shown in equation 1141 



3. Addition of Pulsars 

Adding a new pulsar to an existing array increases the array's overall sensitivity. Equation [T71 
shows how the array's sensitivity increases as a function of a new pulsar's sky location and timing 
residual noise RMS. While we do not have the freedom to choose where we will find the next good 
millisecond pulsar, we do have the freedom to choose where we will look. With this in mind we 
consider how the sensitivity of the NANOGrav PTA would be increased by the addition of a single 
pulsar whose timing residual noise of 200 ns RMS is equal to the current array's median. 

Figure [2] shows the sensitivity Uoveraii as calculated using equation [T7] as a function of the 
location of an addition to the NANOGrav PTA of a single pulsar with 200 ns timing noise RMS. 
As is confirmed by the figure, an additional pulsar will improve the PTA sensitivity regardless of its 
location on the sky, but the improvement represents less than a 6% increase in sensitivity volume 
(the volume of space from which we can detect sources) in all cases. However, some pulsar locations 
will improve it more than others. An additional pulsar in the region in which the PTA is already 
most sensitive improves the sensitivity volume the greatest. This seems appropriate considering 
that the volume of sensitivity goes as d^ . So, if a distance that is already large is doubled, the 
volume increase will be a larger factor than if a small distance were doubled. However, the volume 
sensitivity varied by only 6% as we moved the additional pulsar all over the sky, so it may be wisest 
to search where one is most likely to find pulsars, such as in the galactic plane. One aspect which 
has not been addressed in this manuscript is the coherence of the GW signal between pulsars, i.e. 
the fact that pulsars in similar directions in the sky will show higher correlation between their GW 
signals than those in different parts of the sky (see eq. [T2]) . In addition, in order to confirm that 
the detected signal is a GW and not, for example, an error in earth's ephemerides, or a terrestrial 
clock, we will need pulsars in different parts of the sky to confirm that the spatial signature of the 
detected signal is quadrupolar in nature and not dipolar (ephemerides error) or monopolar (clock 
error). Further study which includes these considerations is necessary to determine the optimum 
strategy. 




Fig. 2. — Shown, the relative sensitivity improvement i^overaii (equation I17p of the NANOGrav PTA 
caused by the addition of a pulsar, as a function of the location of that added pulsar, in equatorial 
coordinates. Every pixel represents the relative sensitivity improvement caused by an additional 
pulsar at that location in the sky. 

4. Optimization With Time Constraints 

The NANOGrav PTA at the time of writing involves 19 pulsars whose timing residuals range 
from 54 ns to 2.2 fis. Monitoring each pulsar requires some fraction of the available observing time, 
which is a valuable resource. How should the available observing time be distributed among the 
different pulsars to optimize the NANOGrav PTA's sensitivity? To explore this question we divide 
the NANOGrav PTA pulsars into "low-noise" (timing residual noise less than 200 ns) and "high- 
noise" (timing residual noise greater than 200 ns) groups, and evaluate the array sensitivity when 
we increase the fraction of observing time spent on low-noise pulsars at the expense of high-noise 
pulsars, and vice versa. 

For a typical timing array pulsar the RMS timing residual noise, a, is i nversely propo rtional 



to the square-root of the time spent observing the pulsar, i.e., a cc t^^'"^ (jVerbiestl l2009l ) . The 
current NANOGrav timing program spends an approximately equal amount of time on each pulsar 
in the array. We consider two alternatives: spending twice as much observing time on pulsars 
in the high- noise group as on pulsars in the low- noise group, and the opposite case of twice as 
much observing time on pulsars in the low- noise group as on pulsars in the high- noise group. This 
translates to the observing time of one group being multiplied by 3 while the observing time of 
the other group is multiplied by |. This, as RMS and t are related above, results in an decrease 
in RMS of the first group by \/3/4 and an increase in RMS of the second group by \/o/2. (In 
both cases the time spent observing the single pulsar with timing residual noise of exactly 200 ns 
is left fixed.) Figure [3] shows the change in volume sensitivity enacted by these adjustments as a 
function of GW propagation direction, s. What we plot is t'(s) where the reference array in the 
denominator (see equation [T8|) consists of the NANOGrav pulsars shown in Table [TJ The array 



used in the numerator is one in which the adjustments described above have been assumed. 

As is noticeable from Figure [3] (with two scales required for plot structure), bettering the 
already good pulsars improves the overall array sensitivity volume by a factor of 1.5 {^overall = 1-5, 
see equation [T71 for definition of Voveraii), while bettering the bad pulsars has quite the opposite 
effect actually worsening the current sensitivity volume {uoveraii = 0.6). The key to understanding 
this result is in noting that SNR /?, residual response r and RMS a are related to each other 
approximately as follows 

(21) 



T 

pcx -, 
a 



and T = TQ{dQ/d) where tq is the amplitude of the residual at earth when the source is at distance 
do- Therefore 

(22) 



To 
POC — , 

ad 



or for a fixed SNR p required for detection, the distance out to which we could detect a source, d, 
is inversely proportional to the RMS a 

d oc l/o-. (23) 

The argument then is similar to that which we made in §3, that volume sensitivity goes as d^ so 
halving an already small RMS a increases the volume sensitivity by a much larger factor than 
halving a larger RMS a. 





0.825 0.83 0.835 0.84 0.845 0.85 0.855 0.86 0.865 



1.12 1.125 1.13 1.135 1.14 1.145 1.15 



Fig. 3. — Sensitivity improvement compared to the current NANOGrav pulsars and current 
NANOGrav pulsar RMS timing values as a function of GW propagation direction s. In both 
panels we plot P(s) as defined by equation 1181 and the text immediately thereafter. On the left we 
improve the RMS of the high-noise pulsars and obtain a volume sensitivity, ^overall = 0.6. On the 
right we improve the RMS of the low- noise pulsars and obtain a volume sensitivity, Voveraii = 1.5. 
Diamonds are the low-noise pulsars, circles are the high-noise pulsars. The median pulsar whose 
RMS we did not alter is a pentagram. All coordinates are equatorial. Please note the scales of the 
two diagrams are different to allow the reader to see the structure of the plots. 
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We can use the same formalism to determine the value of adding more pulsars to the array and 
can use it to ascertain the value of the International Pulsar Timing Array (IPTA) a collaboration 
formed of NANOGrav, the PPTA and the the European Pulsar Timing Array (EPTA). In tabled] 
below the NANOGrav pulsars we have shown the Parkes Pulsar Timing Arr ay pulsars, no t inclu ding 



those overlapping with NANOGrav, as published in the status paper bv iHobbs et al.1 (J2009|). If 
solely PSR J0437-4715 is added to the NANOGrav array, the volume sensitivity increases by 7%. If 
all the PPTA pulsars listed are added to the NANOGrav array, the volume sensitivity increases by 
10%. If instead we imagine an improved situation in which the PPTA pulsars listed are added with 
200ns RMS, with the exception of PSR J0437-4715 which we add at its actual RMS of 100 ns, the 
volume sensitivity is increased by 35%, ie 35% more volume of space is sampled for the same GW 
source type. At the time of writing the European Pulsar Timing Array (EPTA) pulsars for which 



RMSs are available ljFerdman et al.ll2010l ) are already included in this list. Their typical RMSs are 
slightly higher than the values shown here, but there are several reasons to expect that within a 
year the EPTA RMSs will be markedly reduced and that the pulsars will contribute significantly 
to this list. First, the Large European Array of Pulsars (LEAP) project which expects first light 
in late 2010 will create a "tied-array" mode for the 5 European 100-m class dishes into a single 
instrument, riv aling the sensitiyity of Arecibo, but with larger sky coverage. Furthermore, the data 



reported on by lFerdman et al.l (|2010l ) represent new instrumentation at all 5 telescopes, so one can 



expect significant improvement with characterization and optimization of those instruments. 



5. Discussion and Conclusion 

Given the goal to directly detect GWs, the optimization of pulsar timing arrays should be 
well understood. Figure [2] suggests that for the sake of GW detection we may not need to "fill in" 
regions of the sky currently devoid of PTA pulsars, but that rather a clustering of good pulsars in 
one region yields the gratest number of detectable GW sources. The figure shows the sensitivity 
(volumetric) gain produced by adding a new pulsar to the array as a function of the location of 
the added pulsar. The current pulsars in the array are clustered around 18h right ascension and 
declination (as shown by green diamonds on the figure) and in fact the greatest improvement in 
the volume to which the array is sensitive is produced by adding a new pulsar near that already 
existing cluster of pulsars, although the range of improvement is modest in all cases (from improving 
it not at all to improving it by 6%). Here we suggest that if the goal is to maximize the volume of 
space to which we are sensitive to GW sources the location with highest concentration of pulsars is 
slightly favored over other locations, but does not make a significant difference. So perhaps efforts 
to "fill in the gaps" in the spatial arrangements of PTAs are unfounded. However, any new pulsar, 
regardless of its location, is beneficial to the scientific community. Continuing research is needed 
on this particular area to optimize pulsar searches for PTA goals. 

Our initial findings indicate that the sensitivity to burst and continuous GWs can be signifi- 
cantly improved by observing longer the pulsars for which we already have good timing values. In 
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other words, it is found that the intuitive thing to do, observing longer the pulsars that we do not 
have good timing values for, significantly decreases our sensitivity to GWs by almost 50% (Figure 
[3]). This suggestion taken to an extreme yields the ridiculous result that it is best to spend all 
observing time on a signal pulsar. The sensitivity plot in this case would be the familiar beam 
pattern of a single pulsar, but if somehow GWs could be detected with a single pulsar, spending all 
our time on this one best pulsar, assuming the RMS reduces as the square-root of observing time, 
would in fact maximize the volume of space to which we are sensitive. We discussed in ^ that 
this strategy neglects issues with confirming the detection as distinct from clock, ephemerides, and 
other errors, so more work must be done to optimize observing strategies, but this work clearly 
indicates that convention may not be the best means to a discovery. In particular, as the PTA 
is being optimized, what should the observing strategy be? Is the answer different depending on 
whether our aim is to detect single sources or a stochastic background. Presented here is a crucial 
first step toward answering these questions. 

This research was carried out at both Franklin and Marshall College and Pennsylvania State 
University. We wish to express gratitude for the use of facilities at each. We also wish to thank 
Matt Kinsey for his MATLAB interface to HEALPix, and Hannah Williams and Diego Menendez 
for their assistance in developing a web interface to the sensitivity calculations for an arbitrary 
PTA. LSF and ANL gratefully acknowledge the support of National Science Foundation grants 
PHY 06-53462 and AST CAREER 07-48580 respectively. 
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NANOGrav Pulsars^ 





Pulsar 


RMS(;Us) 


1 


J0030+0451 


0.300 


2 


J0218+4232 


0.830 


3 


J0613-0200 


0.110 


4 


J1012+5307 


0.540 


5 


J1455-3330 


0.960 


6 


J1600-3053 


0.190 


7 


J1640+2224 


0.110 


8 


J1643-1224 


1.100 


9 


J1713+0747 


0.055 


10 


J1738+0333 


0.200 


11 


J1741+1300 


0.140 


12 


J1744-1134 


0.130 


13 


J1857+0943 


0.066 


14 


J 1909-3744 


0.054 


15 


J1918-0642 


0.960 


16 


J1939+2134 


0.080 


17 


J2019+2425 


0.910 


18 


J2145-0750 


0.750 


19 


J2317+1439 


0.369 


PPTA Pulsars not part of NANOGrav^ 


20 


J0437-4715 


0.10 


21 


J0711-6830 


1.00 


22 


J1022+1001 


0.50 


23 


J1024-0719 


1.00 


24 


J 1045-4509 


1.00 


25 


J1603-7202 


0.50 


26 


J1730-2304 


1.00 


27 


J1732-5049 


1.00 


28 


J 1824-2452 


1.00 


29 


J2124-3358 


1.00 


30 


J2129-5721 


1.00 


31 


J2145-0750 


0.30 



Table 1: The NANOGrav pulsars with the current RMS timing values at the time of writing. ^This 
is a particular characterization of these NANOGrav pulsars based on communications at the time of 
this writing with the North American Nanohertz Observatory for Gravitational Waves. It is not a 
definitive characterization. We are not presenting the data associated with these pulsa rs but rather 



using them as an example of a realistic PTA. ^RMS values for EPTA pulsars from IHobbs et al 



(|2009|). 



